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We investigate some surfaces of a paradigmatic sp bonded metal — namely, Al(llO), Al(lOO), and 
Al(lll) — by means of the "electron localization function" (ELF), implemented in a first-principle 
pseudopotential framework. ELF is a ground-state property which discriminates in a very sharp, 
quantitative, way between different kinds of bonding. ELF shows that in the bulk of Al the electron 
distribution is essentially jelliumlike, while what happens at the surface strongly depends on packing. 
At the least packed surface, Al(llO), ELF indicates a free-atom nature of the electron distribution 
in the outer region. The most packed surface, Al(lll), is instead at the opposite end, and can be 
regarded as a jellium surface weakly perturbed by the presence of the ionic cores. 



I. INTRODUCTION 

For every given solid surface, breaking of periodicity 
in one dimension will result in a change in the electronic 
states near and at the surface, since the lack of near- 
est neighbours on one side of the surface atoms causes a 
sensible local rearrangement of the surface structure and 
chemical bonds. We investigate here aluminum, which as 
a bulk material is a paradigmatic jelliumlike metal: co- 
hesion and bonding are dominated by electron-gas fea- 
tures, where the ions can be considered roughly speaking 
as a perturbation. What happens at an Al surface, 
however, is much less intuitive: does the surface behave 
essentially like a jellium surface, or do surface atoms play 
a preminent role? Here we investigate this issue and we 
show how a sharp answer is provided by a tool which is 
a very innovative one in condensed matter physics. So 
far, bonding features at a crystal surfaces have been in- 
vestigated by studying either the charge density or the 
local density of states. Instead, we are using here the so- 
called "electron localization function" (ELF), which 
has become recently very popular in the field of quantum 
chemistry. 

We study here three basic choices for the orientation 
of the Al surface: (110), (100), and (111), schematically 
shown in Fig. |l|. These high-symmetry surfaces have a 
rather different packing of the surface atoms: this is also 
visible in Fig. |l|, where the same scale has been used for 
the three surfaces. Some correlation between packing and 
bonding properties is obviously expected, but what is sur- 
prising is that the three chosen surfaces span the whole 
range of possibilities, with (110) and (111) being at the 
two very extreme ends: while the Al(llO) surface promi- 
nently displays well characterized Al atoms, Al(lll) is 
essentially a weakly perturbed jellium surface. The atom- 
iclike vs. jelliumlike character of the electron distribution 
is perspicuously shown by the immediate graphical lan- 
guage of ELF. 

Several appealing features make ELF the tool of choice 
in the present study: ELF is a pure ground-state prop- 
erty, as the density is, but it "magnifies" by design the 



bonding features of a given electron distribution. Fur- 
thermore, ELF is dimensionless, and allows to compare 
the nature of bonding on an absolute scale. ELF pro- 
vides in a very simple way a quantitative estimate of the 
"metallicity" of a given bond, or more generally of a given 
valence region of the system. From our ELF calculations 
it is clear how in the bulk metal the valence electrons 
display a free-electron nature: actually, outside the core 
region, the charge density basically behaves as a free- 
electron gas. This jelliumlike behaviour of the crystalline 
system must be contrasted with the opposite extreme of 
the isolated atom, where the valence-shell region is char- 
acterized by very different ELF values. Through a series 
of ELF contour plots, we will show in the following how 
this free-atom behaviour is almost entirely reproduced 
by the topmost atoms of the least packed metal surface 
studied, namely Al(llO). We stress that ELF discrimi- 
nates an atomiclike valence-electron distribution from a 
jelliumlike distribution in a sharp quantitative way. 

We are adopting here a fully ab-initio method, which 
has become the "standard model" in modern first- 
principles studies of simple metals, covalent semiconduc- 
tors, simple ionic solids, and many other disparate ma- 
terials: j|] namely, density functional theory (DFT) |^ 
with norm-conserving pseudopotentials. This frame- 
work is the natural choice for a real solid system as the 
one chosen for this work, although the ELF has been orig- 
inally introduced and mostly applied, in the quantum- 
chemistry litterature, as an all-electron tool. Anyhow, 
we will demonstrate that the use of the pseudopotential 
approximation makes the ELF particularly meaningful, 
since getting rid of the core electrons and focussing solely 
on the bonding electrons notably simplifies the final out- 
come. 

In a recent very interesting paper. Fall et al. |^ have 
investigated the trend in the Al work function for the 
same three orientations as considered here. They provide 
an explanation of the trend in terms of "face-dependent 
filling of the atomiclike p states at the surface". The 
ELF is an orbital-independent tool, yet it provides a 
concomitant message: we show that there is indeed a 
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face-dependent filling of the electronic states localized in 
the surface region. 

The present paper is organized as follows. In Sec. II we 
introduce the ELF definition, reviewing its foundamcntal 
properties and features. In Sec. Ill we discuss the use of 
the pseudopotential approximation for our ELF analysis. 
In Sec. IV we present our results carried out for the 
three aluminum surfaces. Finally, in Sec. V we draw a 
few conclusions. 



II. DEFINITIONS 

The ELF has been originally proposed by Becke and 
Edgecombe, ||] as a convenient measure of the parallel- 
spin electron correlation. Starting from the short-range 
behavior of the parallel-spin pair probability, they de- 
fined a new scalar function, conveniently ranging from 
zero to one, that uniquely identifies regions of space 
where the electrons are well localized, as occurs in bond- 
ing pairs or lone-electron pairs. This tool has immedi- 
ately shown its power to visualize the chemical bonding 
and the electron localization. The ELF is defined as 



ELF 



1 + [D{v)/Dh[v) 



(1) 



in which D{v) and Dh{r) represent the curvature of the 
parallel-spin electron pair density for, respectively, the 
actual system and a homogeneous electron gas with the 
same density as the actual system at point r. By defini- 
tion, ELF is identically one either in any single-electron 
wavefunction or in any two-electron singlet wavefunc- 
tion: in both cases the Pauli principle is ineffective, and 
the ground wavefunction is nodeless or, loosely speaking, 
"bosonic". In a many-electron system ELF is close to 
one in the regions where electrons are paired to form a 
covalent bond, while is small in low-density regions; ELF 
is also close to one where the unpaired lone electron of 
a dangling bond is localized. Furthermore, since in the 
homogeneous electron gas ELF equals 0.5 at any density, 
values of this order in inhomogeneous systems indicate 
regions where the bonding has a metallic character. 

Savin et al. |^ have proposed another illuminating in- 
terpretation of the same quantity, showing how D can be 
simply calculated in terms of the local behaviour of the 
kinetic energy density, thus making no explicit reference 
to the pair distribution function. This is particularly con- 
venient in our case, since the pair density is outside the 
scope of DFT. Owing to the Pauli principle, the ground- 
state kinetic energy density of a system of fermions is no 
smaller than the one of a system of bosons at the same 
density: ELF can be equivalently expressed in terms of 
the extra contribution to the kinetic energy density due 
to the Pauli principle. The Savin et al. reformulation 
also provides a meaningful physical interpretation: where 
ELF is close to its upper bound, electrons are strongly 
paired and the electron distribution has a local "bosonic" 
character. 



Following Ref. D[y) is the Pauli excess kinetic en- 
ergy density, defined as the difference between the kinetic 
energy density and the so-called von Weizsacker kinetic 
energy functional: jl^ 



Div) = -VrVr'p(r,r') 



1 |Vn(r) 



n(r) 



(2) 



where p is the one-body reduced (spin-integrated) den- 
sity matrix. The von Weizsacker functional provides a 
rigorous lower bound for the exact kinetic energy den- 
sity and is ordinarily indicated as the "bosonic" ki- 
netic energy, since it coincides with the ground-state 
kinetic energy density of a non-interacting system of 
bosons at density n(r). Therefore, D is positive semidef- 
inite and provides a direct measure of the local effect of 
the Pauli principle. The other ingredient of Eq. (j^) is 
Dh{r), defined as the kinetic energy density of the ho- 
mogeneous electron gas at a density equal to the local 
density: 



Dhir) 
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(3) 



It is obvious to verify that the ELF is identical to one 
in the ground state of any one- or two-electron systems, 
while is identical to 0.5 in the homogeneous electron gas 
at any density. 

As commonly done in many circumstances — including 
in other ELF investigations |9[| — we approximate the ki- 
netic energy of the interacting electron system with the 
one of the noninteracting Kohn-Sham (KS) one. We 
therefore use the KS density matrix: 



p(r,r') = 2^</.f^(r)</.:^^(r'). 



(4) 



where <j)f^^{r) are the occupied KS orbitals. Such ap- 
proximation is expected to become significantly inaccu- 
rate only in the case of highly correlated materials. 



III. ISOLATED PSEUDO-ATOM 

The original ELF definition is an all-electron one, and 
has the remarkable feature of naturally revealing the en- 
tire shell structure for heavy atoms. Such a feature is of 
no interest here, since only one electronic shell (the sp 
valence one) is involved in the bonding of Al atoms in 
any circumstances. The pseudopotential scheme simpli- 
fies the landscape, since only the electrons of the relevant 
valence shell are dealt with explicitly: the ELF message 
comes out therefore much clearer, with basically no loss of 
information. In the spatial regions occupied by core 
electrons, the pseudo-electronic distribution shows a de- 
pletion and ELF assumes very low values. Outside the 
ionic cores, in the regions relevant to chemical bonding, 
the norm conservation endows the pseudocharge density 
with physical meaning, as widely discussed in the modern 
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pseudopotential literature. The pseudo ELF carries 
therefore — in the material of interest here — the same in- 
formation as the all-electron ELF, while it removes irrel- 
evant and confusing features due to the inner, chemically 
inert, shells. 

For an isolated pseudo-atom we have by construction 
only a single valence shell, clearly displayed by a single 
and very prominent ELF maximum. In Fig. ^ we report 
our result carried out for an isolated Al atom: in the 
picture one immediately remarks the spherical region as- 
sociated to the ELF maximum (0.91), where the charge 
density is almost bosonic. The black cloud in Fig. ^ in- 
dicates the strong localization of the electronic valence 
shell and perspicuously distinguishes the free aluminum 
pseudo-atom from the crystalline one. 

Indeed looking at the following Figs. ^-|^, and inspect- 
ing the bulk regions of them, the ELF plots quantitatively 
demonstrate that the valence electrons have a predomi- 
nant free-electron nature, with a jelliumlike (or Thomas- 
Fermi) ELF value. The ion cores only provide an "exclu- 
sion region" for electrons (white circles in the plots), but 
basically do not alter the jelliumlike nature of the elec- 
tron distribution outside the core radii: in this sense we 
may regard the ion cores as a "weak" perturbation. ||] 
Actually, the maximum value attained by ELF between 
nearest-neighbor atoms in bulk Al is only 0.61 and — 
outside the core radii — there is a widespread grey region, 
where ELF is almost constant and close to 0.5, thus in- 
dicating a jelliumlike electronic system. This peculiar 
behaviour will help us in the following to understand the 
bonding pattern occurring at different Al surfaces. 

IV. ALUMINUM SURFACES 

All the calculations in this work use a state-of-the- 
art set of ingredients: a plane-wave expansion of the 
KS orbitals with a 16 Ry kinetic-energy cut-off, a set 
of Monkhorst-Pack |l^ special points for the Brillouin 
zone integration, with a Gaussian broadening of O.OlRy 
and a norm-conserving pseudopotential in fully non- 
local form. ||l3| The calculations for the (111) surface 
were performed using a 9-1-6 supercell (9 planes of Al, 6 
equivalent planes of vacuum) and 37 k-points in the irre- 
ducible Brillouin zone. The corresponding values used for 
the (100) surface are an 8-1-6 supercell and 46 k-points, 
and for the (110) surface an 8-1-8 supercell and 48 k- 
points. We have preliminarily investigated the effect of 
surface ionic relaxation on ELF, and found it negligible. 
Similar insensitiveness was found in Ref. |^ upon other 
surface electronic properties: we therefore present results 
for the unrelaxed surfaces. As shown above in Fig. |], the 
three surfaces have a different packing: the packing in- 
creases from Al(llO) to Al(lOO) to Al(lll). A measure 
of this packing is the number of nearest neighbors in the 
surface plane, called coordination in the following: this 
number is 2, 4, and 6 for the three surfaces, respectively. 



In Fig. ^ we report two ELF contour plots along 
two non-equivalent planes passing through the topmost 
(twofold coordinated) atom of the Al(llO) surface. Along 
a direction orthogonal to the surface and passing through 
the surface atom, our calculated ELF attains the maxi- 
mum value of 0.86, thus indicating that the Pauli prin- 
ciple has little effect. Comparing Fig. |3| to |[ where the 
ELF maximum has the close value of 0.91, we quantita- 
tively see see how much the surface atom actually behaves 
as a free-atom in the outer direction. A slightly different 
electronic distribution is present at the four-fold coor- 
dinated Al(lOO) surface in Fig. ^. In this case the ELF 
maximum — along the same direction — attains the lower 
value of 0.80 and shows therefore an electron distribution 
with a less pronounced atomiclike character. 

The almost black regions around the surface atoms in 
Figs. and ^ are also clearly visible in the top views of 
Fig. pT where the contour plots are drawn in the out- 
ermost atomic plane. In this plane the maxima occur 
midway between nearest-neighbors atoms and assume 
the values of 0.77, 0.74, and 0.67 in order of increasing 
packing. This decrease of electron pairing shows a trend 
from covalent-like to metallic-like bonding between near- 
est neighboring surface atoms; bonding to the underlying 
bulk atoms is instead metallic in all cases. 

Finally, in Fig. || we have the extreme case of the 
six-fold coordinated Al(lll) surface. The absolute ELF 
maximum is only 0.73 and lies in a low-symmetry point, 
slightly off the outermost atomic plane. Along the direc- 
tion orthogonal to the surface and passing through the 
surface atom, the ELF maximum is only 0.65: this be- 
haviour clearly marks an almost smooth decay from the 
bulk — where the electrons exhibit a prevailing jellium- 
like distribution — to the vacuum region. Also in the top 
view of Fig. ^ there is no evidence of electron pairing be- 
tween surface atoms, with a wide spread region of almost 
uniform ELF, close in value to 0.5. We can therefore 
summarize the result by saying that the Al(lll) surface 
is essentially a jellium surface perturbed by the atomic 
cores, in a pseudopotential sense: [|l| the perturbation 
obviously originates an "exclusion region" (white plots 
within the core radii), but besides this has little effect 
on the electron distribution in the bulk and at the (111) 
surface. 

The values of the maxima in the direction orthogonal 
to the surface and passing through the surface atom in 
the three cases (0.91, 0.80, and 0.65, as reported above) 
can be interpreted as a measure of the occupation of the 
atomiclike states protruding from the surface. We recall 
in fact that ELF is identically equal to 1 for any one- or 
two-electron system: high ELF values indicate a strong 
localization of the wavefunction. Our numbers show an 
analogous trend as the one pointed out by Ref. ^ , where 
it is shown that the filling of atomic p± orbitals, protrud- 
ing from the surface, decreases indeed with increasing 
packing, i.e. going from (110) to (100) to (111). 
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V. CONCLUSIONS 



(a) 



In conclusion, we have shown how the ELF at the sur- 
faces of a paradigmatic sp-bonded metal distinctly re- 
veals the rearrangements in the electron distribution due _ _ 
to the surface formation and to the changes in the lo- 
cal coordination. We have shown that ELF, in a mod- 
ern DFT-pseudopotential framework, is a powerful tool 
to investigate — without using any spectral information — 
the chemical bonding pattern: this makes the method 
appealing in order to deal with much more complicated 

systems. Jfcl 
Finally, to put this work in a wider perspective, we no- W 
tice in the current literature a quest for other innovative ^ft'l 
tools which address localization and bonding, and rfi 

whose relationships to ELF have not been investigated #^ 
yet. e 
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FIG. 1. Schematic three-dimensional view of aluminum 
surfaces cut along different crystal orientations, (a) (110) sur- 
face, (b) (100) surface, and (c) (111) surface. Black balls cor- 
respond to topmost atoms, gray balls to second-layer atoms, 
and (for the (111) surface) white balls to third-layer atoms. 
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FIG. 2. ELF for an isolated pseudo-atom. The grey-scale 
is also shown: dark (clear) regions correspond to large (small) 
ELF values. The same grey-scale is adopted in all subsequent 
figures. 
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FIG. 3. Al(llO): ELF contour plot in two non-equivalent 
planes containing the topmost atom and orthogonal to the 
surface. 




FIG. 4. Al(lOO): ELF contour plot in two non-equivalent 
planes containing the topmost atom and orthogonal to the 




FIG. 5. Al(lll): ELF contour plot in two non- equivalent 
planes containing the topmost atom and orthogonal to the 
surface. 
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FIG. 6. ELF contour plot m the plane containmg the top- 
most Al atoms for the three different orientations. 
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